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Abstract: The effects of fluid viscosity on the kinematics of a small swimmer at 
low Reynolds numbers are investigated in both experiments and in a simple model. 
The swimmer is the nematode Caenorhabditis elegans, which is an undulating round- 
worm approximately 1 mm long. Experiments show that the nematode maintains a 
highly periodic swimming behavior as the fluid viscosity is varied from 1.0 mPa-s to 
12 mPa-s. Surprisingly, the nematode's swimming speed (~ 0.35 mm/s) is nearly in- 
sensitive to the range of fluid viscosities investigated here. However, the nematode's 
beating frequency decreases to an asymptotic value (~ 1.7 Hz) with increasing fluid 
viscosity. A simple model is used to estimate the nematode's Young's modulus and 
tissue viscosity. Both material properties increase with increasing fluid viscosity. It 
is proposed that the increase in Young's modulus may be associated with muscle 
contraction in response to larger mechanical loading while the increase in effective 
tissue viscosity may be associated with the energy necessary to overcome increased 
fluid drag forces. 



I. INTRODUCTION 

Many, if not most, living organisms move in the realm of low Reynolds (-Re) number, typ- 
ically defined as Re = pUL/p, where p and p are the fluid density and viscosity respectively, 
and U and L are the organism's characteristic speed and length scale. Due to their small 
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size, microorganisms such as bacteria, sperm cells, and various kinds of protozoa move in 
the low Re number regime. In such a regime, linear viscous forces dominate over nonlinear 
inertial forces 

flag 

and locomotion must result from non- reciprocal motion in order to 
break time- reversal symmetry [4], ; this is the so-called "scallop theorem" |f| . 

More than fifty years ago, G.I. Taylor beautifully demonstrated that a slender body 
could swim at low Re number by generating traveling waves along its body Q. An attrac- 
tive example of such a low Re number swimmer is found in the nematode (roundworm) 
Caenorhabditis elegans. This small free-living swimmer _fw 1 mm long) is a well established 
model organism used extensively for genetic research 8] due in part to its short life cycle 
and fully sequenced genome 9|. The motility of C. elegans results from the coordinated 
activity of 95 body wall muscle (BWM) cells, which are highly similar in both anatomy and 
molecular makeup to vertebrate skeletal muscles 1CJ, Anatomically, such BWM cells 
are arranged longitudinally in four quadrants along the nematode body, with two rows of 
muscle cells (dorsal and ventral) in each quadrant [12J, ll3[. Experimental observations have 
shown that motility of swimming C. elegans results from the propagation of bending waves 
along the nematode's body length [bll . flsT ]. which is an example of non-reciprocal motion. 
These bending waves consist of alternating phases of dorsal and ventral muscle contractions 
driven by the neuromuscular activity of BWM cells The adaptability observed in the 
locomotory gait in C. elegans has been hypothesized to be encoded in sensory and motor 
systems (mechanosensation) which allow the nematode to respond to its own movement in 
different physical surroundings [jjj]. In turn, C. elegans is known to switch between dis- 
tinct forms of locomotion when transitioning between locomotion on solid substrates (e.g. 
crawling) and swimming in liquid environments 151 ] . 

The motility behavior of nematodes is a strong function of its material properties such 
as body stiffness and muscle tonus 16]. While it is generally accepted that during loco- 



motion the nematode's body tissues obey a viscoelastic reaction [20|, [21|, [22j, quantitative 
data on C. elegans ' material properties such as tissue stiffness and viscosity remain largely 
unexplored. Experimental investigations on locomotive behavior in C. elegans have mainly 
focused on quantifying kinematic properties using metrics such as bending frequency, ampli- 
tude, wavelength and nematode's centroid velocity Only recently have 
the nematode's material properties begun to be probed. For example, static indentation 
measurements on anesthetized nematodes using piezoresistive cantilevers 23] have provided 
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experimental values for the effective Young's modulus (E) of the C. elegans' "shell", which 
is composed of the cuticle, the hypodermis, and longitudinal muscles. Values are found to 
be on the order of 400 MPa and are much closer to stiff rubber than to soft tissues. More 
recently, estimates of both Young's modulus and tissue viscosity of C. elegans have been 
obtained using kinematic experimental data combined with a dynamic model of locomo- 
bile nematode's tissue properties (E 4 kPa) 



tion 



24] . The authors find that the values o: 



are much closer to biological soft-tissues 
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261 ] . Despite recent efforts, however, there is 



still a dearth of knowledge on overall mechanical properties of C. elegans' body and their 
dynamic response to the external mechanical load imposed during locomotion. 

In this paper, we investigate the motility of C. elegans as a function of fluid viscosity 
(i.e. mechanical loading) in both experiments and in a model. In addition, the nematode's 
effective material properties including Young's modulus and tissue viscosity are estimated 
as a function of mechanical loading. For the range of viscosities investigated here (1.0- 
12 mPa-s), experiments show that C. elegans swim in a periodic fashion while generating 
traveling waves of body curvature with amplitudes that decay from head to tail. While the 
spatial form of the nematode's swimming gait (body amplitude and wavelength) is nearly 
insensitive to increasing fluid viscosity, its beating frequency decreases to a quasi steady- 
state value (~ 1.7 Hz) while maintaining a nearly constant swimming speed (~ 0.35 mm/s). 
A model based on force and moment (torque) balance is able to capture the main features 
of the experiments. The model is used to estimate both the Young's modulus (E) and tissue 
viscosity (77) of C. elegans. We find that both E and rj increase with increasing applied 
mechanical loading (/i) suggesting that nematodes respond to their external environment 
through dynamic changes in their tissue properties (i.e. muscle tonus). 

II. EXPERIMENTAL METHODS 

The experimental configuration (Fig. [T^i) consists of an acrylic chamber that is 1.5 mm 
wide and 1 mm deep; the chamber is sealed with a thin (0.13 mm) cover glass. Shallow 
chambers are used to minimize three-dimensional motion of the nematode. The chamber 
is filled with an osmotically balanced aqueous solution of M9 buffer js|, which contains 
approximately 5 to 10 nematodes per loading. The motion of C. elegans is imaged using an 
inverted microscope and a high-speed camera at 125 frames per second. The high acquisition 
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rate ensures the ability to resolve small linear displacements along the nematode's body 
between consecutive frames. Wild-type (N2) strains of C. elegans are grown and maintained 
on agar plates using standard culture methods Q. They are fed with the E. coli strain OP50. 
All experiments are performed on hypochlorite synchronized young adult animals. Strains 
are obtained from the Caenorhabditis elegans Genetic Stock Center. 

Figure QJb) shows a snapshot of a typical wild- type C. elegans swimming in a buffer solu- 
tion. Nematode kinematics are characterized using in-house image analysis codes. They are 
used i) to create continuous two-dimensional (2D) skeletons from segmented body shapes 
and ii) to calculate the nematode's center of mass (centroid) position. The nematode is 
tracked over multiple bending cycles, and the skeletons represent the nematode's body cen- 
terlines. Results show that the nematode moves in a highly periodic fashion, as indicated 
by the sweeping motion of its tail (Fig. [lb). In particular, body displacement amplitudes 
of the head sweeping motion are always larger than tail amplitudes. In the buffer solution 
(/i ~ 1.0 mPa-s), head and tail amplitudes are typically 490 fim and 375 /mi, respectively. 

The nematode's swimming speed (U) is calculated by differentiating the nematode cen- 
troid position over time. The beating frequency / = 2ir /u and period T = 1/ f are calculated 
from the measured nematode's curvature n(s,t) 24J. Here, u is the angular frequency. The 
curvature is defined as n(s,t) = d<f)/ds, where <fi is the angle made by the tangent to the 
x-axis at each point along the nematode's centerline, and s is the arc-length coordinate 
spanning the nematode's head (s = 0) to its tail (s = L). 

Fluids of different viscosities /i are prepared by adding small amounts of high molecular 
weight (MW) polymer to the M9 buffer solution. The polymer of choice is carboxy-methyl 
cellulose (CMC, MW=2.1 x 10 6 ). Fluid viscosity is varied by changing the polymer con- 
centration in the buffer solution from ppm to 2000 ppm. The fluid viscosities range from 
1.0 mPa-s to 12 mPa-s corresponding to the buffer and the 2000 ppm solutions, respectively 
(Fig. [2]). Since CMC is a flexible polymer, it is natural to expect viscoelastic effects to arise 



in the fluid 



27. 



281 ] . However, the salt ions in the buffer solution screen the CMC molecules. 



This screening tends to keep the molecule in the coiled state even at moderate strain rates 
(~ 10 s" 1 ). As a result, viscoelastic and strain-rate dependent viscosity behaviors are min- 
imized. Figure [2] shows the fluid shear-viscosity versus shear-rate. Note that the shear 
viscosities are nearly-constant even for the highest concentration solutions. The shaded area 
in Fig. [2] corresponds to the mean shear-rates observed in the swimming experiments. No 
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appreciable first normal stress-difference (Nl) is observed for all solutions investigated here. 

III. KINEMATICS OF UNDULATORY SWIMMING 

The swimming kinematics of C. elegans are investigated as a function of increasing fluid 
viscosity (/i). Figure [3] shows (i) the evolution of the nematode's skeletons over a representa- 
tive bending cycle and (ii) the spatio-temporal evolution of the nematode's body curvature 
k(s, t) in the buffer (Fig. [3h and b) and in the 1000 ppm solution (Fig. |3b and d). Envelopes 
of skeleton position are aligned with respect to their centroid, and show that the spatial 
form of the nematode's swimming gait is nearly independent of fluid viscosity (Fig. [3^ and 
c). The nematode's wavelength is nearly constant as well, with a value of approximately 
one body length for the viscosity range investigated here. 

Figure G^b) and (d) shows approximately nine (9) bending cycles of body curvature k for 
the buffer solution and 1000 ppm cases, respectively. The curvature values are color-coded 
such that red and blue represent positive and negative values of k in mm -1 , respectively. The 
2/-axis corresponds to the non-dimensional position s/L along the length of the nematode's 
body. The contour plots show the existence of periodic, well-defined diagonally oriented 
lines. These diagonal lines are characteristic of bending waves, which propagate in time 
along the nematode body length. The waves of body curvature travel along the nematode 
and the curvature magnitude decays from head to tail. This behavior is observed both in 
the buffer solution (Fig. [3b) and in the 1000 ppm case (Fig. (3H). 

Next, the nematode's swimming behavior is characterized as a function of fluid viscosity 
or mechanical load. Figure H] shows the average swimming speed U (Fig. Hk), the beating 
frequency / (Fig. Hb), and the body displacement amplitudes A (Fig. Hb) as a function 
of fluid viscosity. A population of up to twenty (n=20) individual nematodes is used to 
experimentally measure the above quantities for each fluid viscosity (/x). Surprisingly, we 
find that the nematode's forward speed is nearly independent fluid viscosity. For all fluids, 
the nematode's average speed is approximately 0.35 ± 0.01 mm/s. Such speed corresponds 
to a Reynolds number ranging from 0.03 to 0.4; for the range of fluid viscosities investigated 
here, the Re is below unity. Thus, the motility of C. elegans can be considered to be viscous 
dominated even in the buffer solution. 

By contrast, the nematode's beating frequency (/) decreases with increasing fluid vis- 
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cosity (/i), as shown in Fig. H](b). The values of / seem to reach an asymptotic value of 
approximately 1.66 Hz as fi is increased to a maximum value of 12 mPa-s. The absolute 
change in frequency, however, is relatively small such that a ten-fold increase in /i yields 
an effective reduction of approximately 20% in /, from 2.02 ± 0.04 Hz to 1.66 ± 0.04 Hz. 
The spatial form of the nematode's swimming gait seems to be unaffected over the range 
of viscosities; body displacement amplitudes remain near 0.25 mm (Fig. |4b). Such findings 
suggest that C. elegans maintains a consistent swimming gait as it senses its fluidic envi- 
ronment isj], at least for the range of fluid viscosities investigated here. The data in Fig. H] 
indicates that C. elegans may respond to the increased mechanical load by adapting its 
swimming kinematics to maintain a constant swimming speed. This adaptability suggests 
that C. elegans may not be power-limited [jjj]. 

IV. SWIMMING MODEL 

In this section, a simple model is developed to capture the dynamics of the swimming 
nematode. This model will be later coupled with the kinematic data in order to estimate 
the nematode's material properties including its Young's modulus E and its effective tissue 
viscosity rj as a function of fluid viscosity. We begin by considering the motion of the 



nematode in the limit of low Re numbers 



The ratio of its body length (L) to diameter 



(D) is L/D ~ 15 such that the nematode can be approximated as a slender body. The 
nematode's motion is described in terms of its center-line y(s,t), where s is the arc-length 
along the nematode's body and t is time. The swimming nematode experiences no net total 
force or torque (moments) such that, in the absence of inertia, the equations of motion 
become 

OF 

— = C t u t + C n u n , (1) 
^ = -[F,cos(0)-F x .sin(0)]. (2) 

In Eq. (CD), F(s,t) is the internal force in the nematode, is the drag coefficient experi- 
enced by the nematode, Ui is the nematode velocity, and the subscripts t and n correspond 
to the tangential and normal directions, respectively. The drag coefficients C t and C n are 



obtained from resistive force theory 29(. Due to the finite confinement of nematodes be- 



tween parallel walls, corrections for wall effects on the resistive coefficients are estimated for 
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slender cylinders jl[ [30I . 

In Eq. ([2]), M = M p + M a , where M p is a passive moment and M a is an active moment 
generated by the muscles of the nematode; the active and passive moments are parts of a 



total internal moment 



20 



2 ll ]. The passive moment is described using the viscoelastic Voigt 



model 22j such that M p = EIn + r] p I(dK/dt), where / is the second moment of inertia of the 
nematode cross section and n(s,t) is the curvature along the nematode. Qualitatively, the 
elastic part of the Voigt model is represented by a spring of stiffness E while the dissipative 
part of the Voigt model is represented by a dashpot filled with a fluid of viscosity r] p . More 
detail on the model may be found elsewhere 24]. 

In the present model, it is assumed that two homogeneous effective material properties 
describe the entire nematode body, namely (i) a constant Young's modulus E and (ii) a 
constant tissue viscosity t]. In reality, the anatomy of the nematode is more complex; its 
body plan consists of an outer tube (or "shell") separated from an inner tube by a fluid- 
filled pseudocoelom 23]. In particular, C. elegans has a shell-type hydrostatic skeleton (m)]. 
The outer tube is comprised of various organs including the cuticle, hypodermis, excretory 
system, neurons and BWM cells, while the pharynx, intestine, and gonad form the inner 
tube which remains under pressure. For simplicity, however, the complex internal structure 
of C. elegans is treated as lumped into two single material properties (E and rj). Also, the 
nematodes have tapered body shapes at their head and tail, such that the bending modulus 



(Kb = I\/E 2 + u 2 r] 2 ) varies along the body. In the present model, however, it is assumed 
that the nematode's cross-section is constant such that values of Kb remain constant as a 
function of the nematode's arc length s. 

The active moment generated by the muscle is defined as M a = —(EIn a + r) a Idn/dt), 



where position and time dependent pre 



the nematode and r\ a is a positive constant 32|. A simple form for n a can be obtained by 



erred curvature produced by the muscles of 



assuming that n a is a sinusoidal function of time with an amplitude that decreases from 
the nematode's head to its tail 24j. Note that if rj = rj p — r) a > 0, there is net dissipation 
of energy in the tissue; conversely, if r\ = r] p — r\ a < 0, there is net generation of energy in 
the tissue. For live nematodes, it is expected that rj = t] p — rj a < because the net energy 
produced in the (muscle) tissue is needed to overcome the drag from the surrounding fluid. 

Equations (JT]) and (J2J) are simplified by noting that the nematode moves primarily in the 
x-direction (see Fig. [lb). While C. elegans may generate relatively large body displacement 
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amplitudes (A/L ~ 0.25), deflections of its centerline from the x-axis are assumed to be small 
such that the influence of non-linearities is neglected. The small deflection or displacement 
assumption is necessary in order to obtain a closed- form solution for the curvature n(s,t). 
Assuming small deflections from the centerline, s ~ x and cos(0) ~ 1. This results in a 
linearized set of equations given by 

Differentiating Eq. (jlj) with respect to x and combining with Eq. ((3]), we obtain 

d 2 M dy 

-dx^ + Cn m =0 - (5) 

The boundary conditions are such that both the force and moment at the nematode's 
head and tail are equal to zero. That is, F y (0,t) = F y (L,t) = and M(0,i) = M(L,t) = 0. 
Note that the zero moment boundary conditions at the head and tail imply that EIk(0, t) + 
r)I(dK(0, t) /dt) = ElK a (0, t) and EIk(L, t) + T]I{dn(L, t) /dt) = ElK a (L, t). 

Experiments show that the curvature k has non-zero amplitudes both at the head (x = 0) 
and the tail (x = L). In order to capture this observation, we assume that K a (x,t) is a 
sinusoidal wave with decreasing amplitude of the form K a (x,t) = Qo cos ut + Qix cos(ujt — B), 
where Qo, Qi and B are inferred from the experiments 24j. If the curvature amplitude at 
the head is larger than that at the tail, then the nematode swims forward. Conversely, if the 
curvature amplitude is smaller at the head than at the tail, the nematode swims backward; 
if the amplitudes are equal at the head and tail then it remains stationary. 

Equation ([5]) is solved for y(x,t) in order to calculate the nematode's curvature, which 
is defined as n(x,t) = d 2 y/dx 2 . The solution for y(x,t) is a superposition of four traveling 
waves of the general form Ai exp(— (3x cos Pi) cos(f3x sin Pj — ut — <pi) where j3 and Pj can be 
written in terms of id, E, and t], while A t and (pi are constants to be determined by enforcing 
the boundary conditions discussed above (for more details, please see 24J). Note that the 
formulation presented here does not assume a wave functional form for n(x,t). Rather, the 
wave is obtained as part of the solution. 

The nematode body position y(x,t) predicted by the model is fitted to the experimental 
data in order to compute the curvature k. In short, time profiles (about 8 seconds long) of the 
motion are fitted at several points along the nematode body using a Levenberg-Marquardt 
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algorithm. The fitting method yields estimates of the bending modulus (Kb) and the phase 
angle if) = tan -1 (rju/E). The estimated values of K b and if) are later used to calculate the 
nematode's Young's modulus and tissue viscosity. Curvature k contour plots for the buffer 
solution case obtained from experiments and from the model are shown in Fig. 13(a) and 
Fig. [3(b), respectively. For the particular nematode shown in Fig. 5, the estimated values 
of Kf, and if) are 1.87 x 10~ 15 Nm 2 and —78.9°, respectively. The overall estimated values 
of the bending modulus K b and phase angle if), for the buffer solution case, using up to 20 
nematodes is K b = 2.61 x 10~ 15 ±0.22 x 10" 15 Nm 2 and ip = -73.4° ±1.7°. As shown in Fig. 
5(a) and 5(b), the model is able to capture the main qualitative features of the curvature 
profile including the decay from head to tail. 

Next, we compare the temporal evolution of curvature n(t) obtained from the experiments 
and the model as a function of fluid viscosity fi at a single body location s/L = 0.5. Values of 
n(t) for the buffer solution, 1000 ppm, and 2000 ppm cases are shown in Fig. [51(c), 5(d), and 
5(e), respectively. The model seems to be able to capture the experimentally measured values 
of the frequency and amplitude of n(t). The relative error between the values of n(t) obtained 
from the experiments and the model is computed as e = 1 — \ ^ mo dei(U)\) / J2i \ K exp(ti)W- 
The values of e are 3.27%, 1.96%, and 2.31% for the buffer, 1000 ppm and 2000 ppm cases, 
respectively. These values indicate that the analytical solution presented here, while not 
perfect, is able to capture the main features observed in experiments. 



V. C. ELEGANS' MATERIAL PROPERTIES 



Next, the effective material properties of nematodes (E and 77) are estimated from the 
definition of the bending modulus and phase angle if). Here, the nematode is modeled as 



an idealized cylindrical geometry 23 



where the effects of internal hydrostatic pressure 



are neglected; this pressure contribution has been recently shown to have only modest effects 
on nematode body stiffness 



231 ]. In particular, we consider the second moment of inertia / 
of the nematode cross section as represented by the C. elegans 1 "shell" which comprises of 
the cuticle, the hypodermis, and longitudinal muscles. This shell-like structure is defined 
as / = 7r(r^ — rf)/4, where the outer radius of the shell corresponds to the mean nematode 
radius r Q pa 35 /im and the inner radius of the shell is assumed to be pa r Q /2 from electron 
micrographs of the nematode cross-section 23] . 
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Results show that, for the nematodes swimming in buffer solution, E = 0.62 ± 0.05 kPa 
and r\ = —177.1 ± 15.2 Pa-s. The estimated value of E lie within the wide range of values 
of tissue elasticity measured for soft tissues, including isolated brain cells (0.1 — 1 kPa) and 
muscle cells (8 — 17 kPa) 25j. The values of i] for live C. elegans are negative because the 
organism's tissues are generating rather than dissipating energy 24, [32]. We note, however, 
that the absolute values of tissue viscosity \r]\ are within the range (10 2 — 10 4 Pa-s) measured 



for living cells 



26 



34]. 



The values of E and r\ for the nematode swimming in fluids of different viscosities /i 
are shown on a log-log phase diagram in Fig. [6l The estimated magnitudes of both E 
and r) increase monotonically over an order of magnitude with increasing fluid viscosity or 
mechanical load; each data point corresponds to at least twenty experiments. The data 
shows an increase in the nematode's muscle tonus as the nematode senses its environment 
and responds with dynamic changes in its material properties. The increase in the values 
of the Young's modulus is associated with the shortening of sarcomeres and higher density 
of muscle cells typical in muscle contraction 35] . The increase in the absolute values of 
tissue viscosity \r]\ is associated with the larger amount of energy necessary to overcome the 
increased fluid drag or the energy dissipation due to higher fluid viscosity fi. 



VI. CONCLUSIONS 



In summary, we have characterized the effects of mechanical loading on the swimming 
behavior of C. elegans at low Re number by varying the fluid viscosity (/i) . Results show that 
over the range of viscosities tested, nematodes swim consistently with a periodic swimming 
behavior illustrating travelings waves that decay from the nematode's head to tail. The 
spatial form of C. elegans' swimming gait is nearly unchanged while its beating frequency 
decreases only slightly with increasing viscosity. Given such kinematics, nematodes are 
nevertheless able to maintain constant their average forward swimming speed. Such findings 
support altogether C. elegans' locomotive adaptability to its surrounding environment. 

By coupling experiments with a linearized model based on force and torque balances, 
we are able to estimate, non-invasively, the nematode's tissue material properties such as 
Young's modulus (E) and viscosity (77). We find that C. elegans behaves effectively as a 
viscoelastic material with E = 0.62 ± 0.05 kPa and rj = —177.1 ± 15.2 Pa-s for the buffer 
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solution case. Both the magnitude of the Young's modulus and the tissue viscosity increase 
over an order of magnitude as fluid viscosity is increased; values remain within the range 
of that measured for living cells. We speculate that the increase in Young's modulus may 
be associated with the shortening of sarcomeres and higher density of muscle cells typical 
in muscle contraction; usually muscle stiffens as larger mechanical loading is applied. In 
addition, we suspect that the increase in the absolute value of tissue viscosity |?7| is associated 
with the larger amount of energy necessary to overcome the increased fluid drag or the energy 
dissipation due to larger fluid viscosity //. 
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Figure Legends 

Figure CD 

(a) Schematic diagram of the experimental setup. The nematode Caenorhabditis elegans 
(C. elegans) is placed in a shallow channel that is 1.5 mm wide and 500 /im deep containing 
an aqueous solution. The channel is sealed with a thin (0.13 mm) cover glass. C. elegans is 
imaged using standard microscopy and a high-speed camera, (b) Visualization of instanta- 
neous body centerline or nematode's skeleton with resulting centroid and tail-tip trajectories 
over multiple body bending cycles. Image acquisition rate is 125 frames per seconds. The 
nematode's undulations lie within the microscope focal plane. 

Figure [2} 

Fluid characterization using a stress-controlled rheometer. Shear-viscosity /i as a function 
of shear-rate for both buffer and polymeric solutions. Shaded area in the plot corresponds 
to shear rates (7) associated with C. elegans swimming motion. For the range of shear 
rates considered, the values of fluid viscosity are nearly constant and fluid solutions may be 
regarded as effectively Newtonian. No appreciable normal first normal stress difference (Nl) 
was observed even for the highest polymer concentration (2000 ppm). 

Figure [2 

Skeletons and spatio-temporal evolution of C. elegans' curvature swimming in different fluid 
viscosities. Examples of the evolution of nematode skeletons over one representative bending 
cycle in (a) the buffer M9 solution (/j, pa 0.8 mPa-s) with T pa 0.49 s, and (c) in the 1000 ppm 
CMC solution (/i pa 3.59 mPa-s) with T pa 0.62 s. The data reveals a well-defined envelope 
of elongated body shapes with a wavelength corresponding approximately to the nematode's 
body length. Envelope shapes are aligned by their centroids and are nearly invariant with 
increasing viscosity, (b) and (d) Representative contour plots of the experimentally measured 
curvature (k(s, t)) along the nematode's body centerline for (a) the buffer solution and (d) 
1000 ppm CMC solution, respectively. Approximately 9 bending cycles are shown over 4 
seconds. Red and blue colors correspond to positive and negative k values, respectively. The 
units of k are in mm -1 . The y-axis corresponds to the dimensionless position s/L along the 
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C. elegans 1 body length where s = is the head and s = 1 is the tail. 
Figure @} 

Nematode tracking measurements as a function fluid viscosity fi. Nematode's (a) swimming 
speed ([/), (b) beating frequency (/), and (c) body displacement amplitude (A). Both the 
forward speed and displacement amplitudes are nearly constant with increasing fluid vis- 
cosity The nematode's bending frequency decreases to a quasi-steady value with increasing 
fluid viscosity or mechanical load. 

Figure [5t 

(a) Representative contour plot of experimentally measured curvature n(s,t) along the ne- 
matode's body centerline for the buffer solution case, (b) Corresponding contour plot of 
curvature values obtained from the model. The model captures the longitudinal bend- 
ing wave with decaying magnitude, which travels from head to tail. For the nematode 
shown, estimates of the bending modulus and phase angle yield Kb = 1.87 x 10 -15 Nm 2 and 
if) = —78.9°, respectively. Such values correspond to E = 0.32 kPa and \rj\ = 134.5 Pa-s. 
(c)-(e) Comparison between experimental and theoretical curves of n(s,t) at s/L = 0.5 for 
the buffer solution, 1000 ppm, and 2000 ppm cases, respectively. The relative error between 
the model and the experiments is less than 4% (see text). 

Figure [6t 

C. elegans' estimated material properties as a function of fluid viscosity. The estimated 
Young's modulus E and absolute values of tissue viscosity \rj\ increase monotonically with 
increasing fluid viscosity. 
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